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Qh! Abstract: Recently Hertog, Horowitz, and Maeda have argued that cosmic censorship 

(— I I can be genericahy violated in string theory in anti-de Sitter spacetime by considering a 

collapsing bubble of a scalar field whose mass saturates the Breitenlohner-Preedman bound. 



We study this system numerically and find that for various choices of initial data black 
holes form 
is upheld. 



■ holes form rather than naked singularities, implying that in these cases cosmic censorship 
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1. Introduction 

Inasmuch as they mark a breakdown of classical general relativity, spacetime singularities 
pose an important challenge to theories of quantum gravity such as string theory. Singular- 
ities are common occurrences in the study of cosmology and black holes, but their correct 
physical description still seems to be beyond our grasp. Penrose proposed the concept of 
cosmic censorship, which states that singularities which are formed by nonsingular initial 
data are always hidden behind a horizon [||]. 

In the case of cosmology it has proven quite challenging to find controlled examples 
within string theory [^] [^] ^ ||] , while for black holes the existence of the event 
horizon makes it difficult to study the physics of the singularity itself p!o|]|ll||l^|13]|14|. 



To make progress it would be very helpful to find a naked singularity within a well 
understood background of string theory. Finding a naked singularity in an anti-de Sitter 
background would be especially promising, as the nonperturbative tools of the AdS/CFT 
correspondence could be brought to bear on the problem. An interesting proposal for 
obtaining such singularities was put forward in a recent paper by Hertog, Horowitz and 
Maeda (HHM) |l|( see also (l|]). Their argument uses the special properties of scalar fields 



which saturate the Breitenlohner-Freedman bound |17|[18|, which describes how tachyonic 
a field in AdS can be without leading to an instability. If one starts with a large and 
homogeneous bubble of such a scalar field the resulting evolution leads to a singularity in 
finite time. To argue that the singularity is naked, HHM claim that the energetics of the 
system is such that no black holes can form. 

In this paper we present results on the numerical study of such collapsing bubbles. 
For all the cases we have studied, we will see that the evolution is consistent with the 
formation of black holes rather than naked singularities. HHM have actually considered 
two distinct scenarios. For the first, one imposes standard boundary conditions in AdS, and 
the argument revolves around the fact that the initial data has negative mass; black holes 
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on the other hand have positive mass, which might seem to prevent them from forming from 
the assumed initial data. However, there is a loophole: this definition of the mass is not 
conserved, as HHM themselves point out. Numerically we will find that the mass rapidly 
becomes positive, and that a trapped surface subsequently forms outside the bubble. In 
the second scenario HHM impose Dirichlet boundary conditions at some cutoff position 
inside AdS. In this case, the mass is conserved, but one also has to consider the possible 
existence of hairy black holes whose mass is lower than that of the initial data (the no hair 
theorem does not apply in the presence of the cutoff). We will compare the mass of the 
black hole with the mass of the initial data for various choices of parameters and find that 
black holes have sufficiently small mass to form. Similarly, our numerical investigations 
of the collapsing bubble with Dirichlet boundary conditions again leads to the conclusion 
that black holes form rather than naked singularities. Thus, for all the examples we have 
studied cosmic censorship is upheld. However, we should note that we have by no means 
exhausted all possible choices of parameters, including some suggested in [|^], and so it is 
possible that cosmic censorship can be violated in some other regime. 



2. Equations of motion 



We consider the action for a scalar field coupled to gravity in 5 dimensions with a negative 
cosmological constant. 



S = jd''x^(^R + Q-]^[dct>f-V{ct>)). (2.1) 



The potential corresponding to a free scalar which saturates the Breitenlohner-Freedman 
bound is 

V{cp) = -2(t?. (2.2) 



For now we will assume |2.2| , although we will also consider the full potential arising in 
supergravity in section 3. We have chosen our units such that the background AdS solution 
is 

ds^ = -{l+r'^)dt^ + -^—^+r'^dn^^, (/)=0. (2.3) 

Spherical symmetric solutions will be considered throughout. The general such metric 
in the ADM decomposition is 

ds^ = -N\t, rfdt^ + L{t, rf{dr + N''{t, r)dtf + R{t, rfdfll (2.4) 

It is very convenient to choose R{t, r) = r so that the coordinate r has a clear physical 
meaning. In our numerical analysis we will be looking for the appearance of trapped 
surfaces, and we now state the condition for their existence. Consider radial and future 
directed null geodesies, obeying 
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A trapped surface exists whenever ^ < for either choice of sign. This can be rewritten 
in a more useful form as follows. Define 

^ = ^' + ^'-l2 + (^J • (2-6) 

\i is proportional to the quasilocal mass, and essentially measures the total amount of 
energy contained within an 5^ of radius r. If we solve (^) for L and substitute into ( p^ ) 
we find 

^n=(-A.'±^(ArO^+(l+.^-ii)(A..)^)^. (2.7) 

We conclude that trapped surfaces occur for 1 + - 4 < 0. For given constant /i, this is 
the same as saying that r is inside the horizon of the AdS-Schwarzschild black hole. We 
will refer to 1 + — < as the "horizon function" . 

{i was defined above in any coordinate system of the form ( |2.4| ) with R = r. In our 
numerical study we will choose to set A^^ = and write the metric as 

ds' = -n{t, rf{l + r'- ^)dt' + ^^^^ + r'dnl (2.8) 



r 



1 + - 



Setting n = 1 and fi = constant gives the standard form of the AdS-Schwarzschild black 
hole, the mass being M = Svr^/x. Since the definition of fi in (2^) agrees with that in (^.6|), 
we again see that trapped surfaces occur for l+r'^ — yi < 0. But given the form of the metric 
(|2.8| ), it is clear that the coordinates will break down precisely when trapped surfaces are 
encountered. Therefore in our numerics we will look for 1 + — ^ becoming very small, 
and assume that the subsequent evolution is sufficiently smooth that trapped surfaces do 
indeed occur. We will also give evidence that the trapped surfaces are associated with black 
hole formation, although it would be interesting to check this directly by using coordinates 
which can extend into the region containing the trapped surfaces. 

With spherical symmetry, the gravitational field has no dynamical degrees of freedom, 
and the metric functions n and /i on a given time slice can be determined in terms of <j) and 
(j) evaluated on the same time slice. In particular, two linear combinations of Einstein's 
equations yield the constraints 



' ^r30^ + ^(l + r^-^)(5+rW), (2.9) 

2 



Here, vr^ is the canonical momentum 



„ 3.n.3+'W)- (2.10) 



TTs = JT^^. (2.11) 

We also need the scalar field equation of motion, which is 

= n(r^{l + r^- ^)0" + (Sr^ + 5/ - /x + ^rV)*/"' + 4r^(/') (2.12) 
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where and (|2.10[) were used. 

For numerical purposes it is useful to redefine variables. We map r £ [0, oo] to x G 
[-1,1] by 

2 _ 1 + X 



1 — X 

The asymptotic behavior of the scalar field can be simplified by defining 



(2.13) 



iP = {1 + r'^)(j) = -^(j) (2.14) 
1 — X 

vr^ = (1 - x)~^/\l + x)-^/\^. (2.15) 

We now discuss the boundary conditions. In solving (|2.9D and ( |2.10D we can choose 
the values of /i and n at the origin, and we choose these to be 

^(t, r = 0) = fi{t, X = -1) = 0, n{t, r = 0) = n{t, x = -1) = 1. (2.16) 

Note then that our time coordinate measures the proper time at the origin. As for the 
scalar field at the origin, we need only demand that ip is nonsingular there. We also need 
a boundary condition on the scalar field at infinity. Solving the free scalar field equation 
in AdS yields the asymptotic behavior 

ct(t) In T ~ 

't'^^+Pi^)— ^ ^{x)^a{t)+(3{t)\n{l-x). (2.17) 

The two asymptotic behaviors correspond to the normalizable and non-normalizable modes 
familiar in the study of fields in AdS. We are allowed to freely specify f5{t) as well as the 
values of (j) and (j) on a-n initial time slice, and then a{t) is determined by the equations of 
motion (and analogously for V', a and j3). Fluctuations around the ordinary AdS vacuum 
correspond to setting f3{t) = (3{t) = 0; on the CFT side of the AdS/CFT correspondence 
this corresponds to evolution with respect to the unperturbed CFT Hamiltonian. We will 
refer to these boundary conditions as the "standard" ones. Besides the standard boundary 
conditions, HHM consider a second set of Dirichlet boundary conditions, which we'll return 
to in the next section. So, expressed in terms of ip, the standard boundary conditions just 
say that should be finite as x ^ 1, but is otherwise unrestricted. 
The equations of motion can now be written 

^ = 2(^2 - \ l^)mr^, (2.18) 

= ^(l-x)(2(l + x)-(l-x)V)nV'"- (l + ^\^^/^)"V' 

+ n(il- 3x) + + + 1 (1 _ ^)2(i + ^)^2\^/ 

V 1 + X / 

- ^(l-x)(l+x)n^^ (2.19) 

where the metric function n can be expressed as follows 

^ = eArfi/(i) (2.20) 
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and we defined 



and is given by 

+-(1 + xf{l - x){ij'f - -(1 + x) VV'je-^-i'^^^^*) (2.21) 
6 3 J 

/ = i ((1 - xf^l + 1(1 - x'm - xW - i^f) . (2.22) 

Note that the only potentially singular terms are of the form but /x ~ (1 + x)^ as 
X — > —1, so this is not a problem. 

Given the way in which ^ appears in the metric, its value at infinity gives a natural 
definition (up to a factor of Svr^) of the total mass associated with the spacetime. However, 
for scalar fields which saturate the Breitenlohner-Freedman bound, it turn out that this 
definition of mass is not conserved. From the definition ( 2.21| ) and the equations of motion, 
one can show 

|Mi,^ = l) = -^|[V'(t,x = l)2]. (2.23) 
A conserved mass can therefore be defined as 

2 

m = n{t, X = 1) + -ij{t, X = 1)2. (2.24) 

Another important fact is that fi can be negative, while m obeys a positivity theorem. 

Following HHM, our initial data will consist of a large homogeneous bubble of scalar 
field, with vanishing time derivative, and with the standard AdS boundary condition at 
infinity. So at t = we have 

2 

tp = (po for X < xq (2.25) 

1 — X 

-0 ~ finite for x ^ 1 (2.26) 
vr^ = (2.27) 

with (pQ being the constant value of (j) inside the bubble. We take the bubble to be suf- 
ficiently large that a singularity is guaranteed to form at the center of the bubble within 
finite time.^ Furthermore, we restrict the size of the bubble and the scalar field profile 
outside the bubble to be such that the asymptotic value of ^ is negative. Were /i conserved 
this initial data would certainly evolve to a naked singularity, since black holes have > 0. 
However, we know that /i is not conserved, and the question becomes whether it can in- 
crease sufficiently rapidly to allow a black hole to form before a naked singularity appears; 
we will see that this is indeed what happens. 

3. Results 



We solve the equations ( p. 18 ) and ( 2.19| ) numerically using C code by evaluating the inte- 



grals for n and fi on a given time slice, and then evolving ip and vr^ forward by a simple 
^See mM for a discussion of the required size. 
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time stepping. We take the number of spatial lattice points to be a few thousand. Stability 
requires that the time step is approximately the square of the spatial lattice spacing, and 
we have taken dt = dx'^/2. We have done the numerical evolution using various number 
of lattice points, and a refinement of the lattice does not show any sign of a numerical 
instability. In the course of the evolution n(x, t) becomes very large near the boundary. In 
order to avoid a numerical instability we have to adjust the time stepping dt — > dt/n{t,x). 
We evolve in time until the numerics break down, which occurs when the horizon function 
1 + r'^ — becomes very small, and this is associated with the appearance of trapped 
surfaces. 



3.1 Standard boundary conditions 



With 4001 lattice points, we have chosen various values for (pQ and xq in the initial data. 
Outside the bubble, the initial data for •i/^ is taken to be linear in x, with slope matched to 
the first derivative of tp just inside the bubble. In the following we will present the results 
for two cases. In the first we present a fairly generic example using the potential and 
in the second we instead use the full supergravity potential. 



Case 1. Quadratic potential: For xq = .583 and (pQ = 0.2, we find that our numerics 
break down at t final = 0.8837. The initial profile of ip{x) is given by Fig. 1 
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Figure 1: Initial profile of i}{x) &i t — 
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whereas the final profile is given by Fig. 2 
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Figure 2: Final profile of ^(x) at < = t final 



The plot of the conserved mass m{t) versus nonconserved mass fi{t) in Fig. 3 reveals that 
indeed m{t) is conserved all the way to tfinai- The fact that m(t) is conserved to great 
accuracy all the way to t = tfinai constitutes a good check of our numerics. 
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Figure 3: Conserved mass versus nonconserved mass as functions of time 
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As discussed in section 2, the vanishing of the horizon function 1 + — ^/r^ indicates the 
occurrence of a horizon (or more accurately an apparent horizon) at Xh = —0.376, as seen 
in Fig. 4 
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Figure 4: Horizon function at t = 


~ ifinal 



The Ricciscalar remains well behaved near x = Xh, although it grows in the region x < Xh- 
At t = t final the Ricciscalar is shown in Fig. 5 
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Figure 5: Ricciscalar at t = t final 
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Case 2. Supergravity potential: The quadratic potential (|2.2| ) for the scalar field (f) is 
an approximation to the exact supergravity potential 



y((/>) = -2exp( 



-^</.)-4exp(--^0)+6 



(3.1) 



Numerically, this is more expensive to evaluate. For the same initial data as shown in Fig.l, 
with xq = 0.583 and = 0.2 and 2001 lattice points, the numerical evolution breaks down 
after t final = 0.8217. The results are very similar to the ones for the quadratic potential. 
The final profile is plotted in Fig. 6. 




Figure 6: Final profile of ipix) at i = t final 

The horizon function at t = TfinaU ^ plotted in Fig. 7, indicates that the horizon forms 
near = —0.41. 
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Figure 7: Horizon function a.t t — t final 

Hence the full nonlinear potential displays the same characteristic features as the quadratic 
potential. 
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These two numerical results are a sample of the generic features of almost all initial 
configurations we have studied. The nonconserved mass rapidly becomes positive and an 
apparent horizon forms. The size of the horizon grows with the conserved mass, and we have 
studied initial data where the horizon radius is considerably larger than 1. We interpret the 
apparent horizon to be associated with black hole formation since it appears well outside 
the bubble in a region where the curvature seems to remain small. In particular, we are 
distinguishing our trapped surface from those trapped surfaces which we know will appear 
in the homogeneous region of the bubble as the singularity is approached; the latter have 
nothing directly to do with black hole formation. 

We also seem to find the first indications for the decay of the scalar field, in line with 
the no-hair theorem, and a likely endpoint of the evolution is an AdS-Schwarzschild black 
hole without scalar hair. However, since our numerics break down when the horizon forms, 
one has to use more sophisticated numerical methods to settle this question. In any case, 
it seems clear that cosmic censorship is not violated in this setup. 

3.2 Dirichlet boundary conditions 

HHM have proposed a second scenario for producing naked singularities, which involves 
imposing Dirichlet boundary conditions at some large but finite radius, r = Ri?' We take 
constant cj) = (f)Q inside a bubble of radius Rq, and set (j){Ri) = (pi- To minimize the energy 
we take (/>(r) to obey the linearized equations of motion outside the bubble: 

(l) = (po {r < Ro) 

^=4 + 6Ml±l!l {Ro<r<Ri) (3.2) 

where a and b are determined by the boundary conditions at Rq and Ri. HHM actually 
consider slightly different initial data which do not satisfy the equations of motion outside 



the bubble, but we consider |3.2| since it will clearly give lower mass and thus be more 
favorable for violating cosmic censorship. 

We will reduce to a three parameter family of initial data by minimizing the mass with 
respect to (po while holding Rq, Ri and (pi fixed. This yields the mass formula 

toi = -l^lRj(^l + Rl{l- (2 + ^2) 1^(1 + ^2) ^ +'rI) ln(l + Rl) + 2Rl ) ) " ^^'^^ 

Following HHM, we further require (/>o ^ 1 so as to be able to neglect backreaction. In 
this regime, HHM show that if a black hole forms then the size of the event horizon on the 
initial surface must be at least (P^^^Rq. 

For Ri ^ Rq ^ 1, which HHM consider, we have /^initial < 0, and fi is conserved due 
to the Dirichlet boundary condition. But to establish a violation of cosmic censorship one 
needs to further argue that any black hole obeying the same Dirichlet condition at and 

2 /3 

with horizon size at least </>q Rq, has /Ubh > /^initial- 



■^We thank Thomas Hertog and Gary Horowitz for discussions regarding the material in this section. 
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We therefore look for static black hole solutions, satisfying 



^J■+^r(p --(1+r -^)r((/>) =0 
r2(l + - 4)'A" + (3r^ + 5r^ - /i + -rV^)0' + 4r3'A = 0. (3.4) 



We take the horizon to be at r = i?^, and so n{Rs) = Rg + R^. The equations of motion 
can be seen to fix (f)'(Rs) in terms of and so the black hole solution is completely 

specified by Rs, Ri and 01. We numerically integrate the equations outward from the 
horizon to and then adjust Rs and (piRs) so as to satisfy = (pi- We also have 

2 /3 

to satisfy Rg > (pQ Rq in order to be sure that our horizon is large enough to enclose the 
region in which we know singularities will form. Although we were able to use the free 
field approximation to compute the mass of the bubble, it turns out that backreaction is 
quantitatively important for the black hole, and so we are forced to solve the full nonlinear 



equations 3.4 



We have carried out this procedure for a few parameter choices, as summarized below. 



Example 1: We take 

Ri = 20000, Rq = 100, (pi = 10"^ (3.5) 



Then (pQ = 0.006, hence the backreaction for the bubble can indeed be neglected. The mass 

2 /3 

of the bubble is then /^initial = —98198.5 and the critical horizon radius is (p^ Rq = 3.429. 
We then succeed in finding a black hole solution with Rg = 4 and /^bh = —100165. Since 
Mbh < Minitiai the energetics allow black hole formation, and so there is no reason to expect 
cosmic censorship to be violated. 



Example 2: We take 

Ri = 100000, Ro = 500, (pi = lO"'^ (3.6) 

2 /3 

yielding /^initial = —613742, and the critical horizon radius is (pQ Rq = 3.699. We then find 
a black hole solution with Rg = d and fi^h = —630601. So we again find fi^h < ^initial- 
It is possible that there exists some other choice of parameters such that ^bh > /^initial; 
but we have not encountered such a case so far. 

Now we turn to the numerical time evolution in the Dirichlet case. Using our variables, 
the regime InRi/lnRo ^ 1 considered by HHM is rather inaccessible due to the large 
hierarchy of scales. We will instead consider similar initial data as we studied in the case 
of standard boundary conditions, and return to the lni?i/lni?o ^ 1 case in a future 
publication. 

For the numerical evolution we chose initial data with xq = 0.583, (po = 0.2, and 
with 3981 lattice points. The Dirichlet boundary condition is imposed at xi = 0.99. The 
numerical evolution breaks down at t final = 0.9139. 
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The initial profile of iIj{x) is shown in Fig. 8 
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Figure 8: Initial profile of tp{x) at t = 



and the final profile is shown in Fig. 9. 
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Figure 9: Final profile of ^(x) at t = t final 
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The horizon function 



at t = t final is plotted in Fig. 10, and the Ricciscalar is weh be- 



1.4; 

1.2- 



1 




-0.9 -0.8 -0.7 -0.6 -0.5 -0.4 
Figure 10: Horizon function a.t t ~ t final 



haved near the horizon. We have also checked that the mass is conserved to excellent 
accuracy. The approximate vanishing of the horizon function indicates that the evolution 
with Dirichlet boundary conditions leads to a singularity shielded by a horizon instead of 
a naked singularity. 



4. Discussion 



We have given evidence that cosmic censorship survives the challenge from collapsing bub- 
bles of Breitenlohner-Freedman mass scalar fields in AdS. In particular, the existence of 
negative mass for the initial data is consistent with black hole formation since the mass 
becomes positive rather early in the evolution. In the case of Dirichlet boundary condi- 
tions there is still a possibility of violating cosmic censorship since we have not exhausted 
the full parameter space. However, we should point out that if one truly wishes to form a 
naked singularity in string theory then we are presumably restricted to the standard bound- 
ary conditions; being a quantum theory of gravity, string theory only admits asymptotic 
boundary conditions and not Dirichlet conditions at finite locations. 

There have been a number of previous numerical studies of the dynamics of scalar 
fields coupled to gravity in AdS spacetimes |19|[20|. In these papers the main focus was 
on massless scalar fields. The study of gravitational collapse and black hole/singularity 
formation using numerical techniques is a very interesting problem in general. 

In the future we hope to study various different scalar field potentials, for example 
those considered in |l^ |21[ | , and to employ coordinates which continue through the horizon 
of the black hole. Our numerics were based on the simplest discretization and evolution 
schemes, and it would be interesting to apply more sophisticated techniques like mesh 
refinement or excision to this problem and check our results. In particular it would be 
interesting to investigate the critical phenomena associated with gravitational collapse (see 
1 22] 1 23] for reviews) in the context of the AdS/CFT correspondence. 
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